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Abstract 

We study the production asymmetry of charm versus anticharm mesons in photon- 
proton interactions. We consider photon gluon fusion plus higher order corrections in 
which light quarks through vector meson -proton interactions contribute to the cross 
section. Non perturbative effects are included in terms of a recombination mechanism 
which gives rise to a production asymmetry. 



This work was supported by: CLAF (Brazil) and CONACyT (Mexico), CIEA-JIRA (Mexico) 



1 



1 Introduction 



In high energy photon-hadron interactions, charm production is expected to be dominantly 
produced by photon-gluon fusion processes. According with QCD perturbative calculations, 
this mechanism produces in equal amounts charm and anti-charm. 

However, recent measurements of charm meson production [0 indicate that there are impor- 
tant nonperturbative QCD phenomena in the production process that induce an asymmetry 
in charm and anti-charm production. 

This phenomena has been observed in hadron hadron collisions [|1 and is well known as 
"leading effect". It has been the subject of many models of particle production and several 
mechanisms have been proposed to explain it ||, Here we study the xp distribution 
of D"^ and mesons produced in photon-proton collisions in the framework of a two- 
components model that has been used before to successfully describe the asymmetry in pion 
proton interactions Q]. 

The production of D mesons in the model is assumed to take place via two different 
processes, namely QCD parton fusion with the subsequent fragmentation of quarks in the 
final state and conventional recombination of valence and sea quarks present in a vector 
meson fluctuation of the photon. 

The asymmetries obtained with the conventional soft charm component as well as with a hard 
charm component in the photon, are presented. We compare both with the experimental 
data available. 

To quantify the difference in the production of charm and anti-charm mesons an asym- 
metry A is defined as in 

where Nc and iVg are the production yields. The asymmetry has been observed to be a 
function of both xp and the transverse momentum pt. 

This paper is organized as follows: the photon gluon mechanism for charm production 
will be discussed in section 2. This mechanism is not responsible of a production asymmetry 
between charm and anti-charm. In section 3 we discuss the contribution of the resolved 
photon in the frame of a Vector Dominance Model (VDM) component. In subsection 3.1 we 
calculate the QCD interaction of the resolved photon. In subsection 3.2, a recombination 
mechanism is discussed. While the QCD production mechanisms in photoproduction are 
the same for charm and anti-charm production, recombination favors the formation of 
over and over D^. In section 4 the various components are put together and the 
asymmetry is estimated. Finally some conclusions are drawn in section 5. 
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2 The photon-gluon fusion mechanism 

In this section we outhne the calculation of the photon-gluon fusion contribution at Leading 
Order (LO) in pQCD to the D- meson inclusive xp distribution. 

The processes involved in the photoproduction of charm at LO in pQCD are depicted in 
Fig. 1. The corresponding formula in the parton model is 

where g{x, Q^) is the gluon probability density in the proton and 

^^^^"rfW^ = 2i^(2;0^('") ^^P- + p. - - p.) \M\ . (3) 

The squared invariant matrix \Mf in terms of the Mandelstam variables is given by (see 
e.g. B i) 



1 / mg -t ^ ml-u \ ^2 / ml ^ ml \ ml ^ ml V 

2\ml — u ml — ij \ml — i ml — u J Img — t ml — u J 



(4) 



and 



s = xs 

i = ml — mrVse'^" 

u = ml — xmT^/se^" , (5) 

where s is the cm energy of the 7 — p system, m"^ = ml x is the momentum fraction of 
the gluon inside the proton and i/c = (1/2) ln[{Ec — Pc) / {Ec + Pc)] the rapidity of the c-quark. 

After integrating out Eq. ^ on the c-quark variables and the momentum fraction x, the 
inclusive cross section for the production of charm (anti-charm) is given by 

da _ \/i r xgi.x.Q'^) da 



with 



dxp 2 J ^ E{y/s - mrey-) di 



da iraeasiQ^) .^.2 



dt s 
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In Eq. 1^, as is given by 



X = . (7) 



(33 - 2Nf)log^ 



with Nf = 4 and = A4 according to the gluon distribution used in Eq. |[ 

In our calculations we use the GRV-LO gluon distribution with = 4m^, rric = 1.5 
GeV. The D- meson inclusive xp distribution including fragmentation is given by 

dan f dz da etc) „ , . 

- ' '-Dj:)/^{z) 



X 

We use the Peterson fragmentation function, 



dxp J z dx 



Dd/c{z) = ^. (10) 

z 1 - i - 

[ 2 1 — Z 

Inclusion of Next to Leading Order (NLO) contributions into the D-meson cross section 
does not produce appreciable changes neither in the form of the D-meson distribution or in 
its normalization (see e.g. At NLO the cross section for the production of an anti-quark 
differs from the cross section for the production of a quark, but this effect is small @]. 

3 Light quark corrections to charm photoproduction 

We will identify the hadron like part of the photon with p and u vector mesons neglecting 
the contribution of heavier resonances which have smaller couplings to the photon. The p 

and the u can be regarded as two states systems 

— ' ■ ■ ' ^ : yaa — aa j . 

(12) 



p = —1= [uu + dd^ ; 1^ = yi-"^ ~ ddj . (11) 



The photoproduction of D mesons may take place by QCD interaction of the vector meson V 
with the proton and/or by recombination of its constituent quarks (see Fig. 2). Therefore, 
in obtaining the differential cross section of the process Vp DX, we will consider two 
possible processes, 

^^.(Kp ^ DX) = '^^^{Vp ^ DX) + ^^^{Vp ^ DX). (13) 



dxp dxp dx} 

In order to calculate these two contributions to the total cross section, we assume that 
the momenta distribution of the quarks of a p, u; meson are the same than inside a pion. We 
will use the GRV parametrization for the parton distribution in the pion. 
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3.1 QCD resolved photon contribution 

The photon may interact through ist constituents with the partons in the proton. In the 
parton fusion mechanism, D^{D^) mesons could be produced via the qq{gg) — > cc with 
the subsequent fragmentation of the c(c) quark. The contribution of the hadronic (or re- 
solved) component of the photon is given by the usual formula that describes hadron-hadron 
interactions, 

where Xi/j(xi, /i^) is the parton distribution in the resolved photon, X2fi{x2, /i^) is the parton 
distribution in the proton, E is the energy of the fragmenting charm quark and Do/dz) is 
the fragmentation function. 

The partonic cross section in Eq. |14| is given by 

da 7ra^(yu) 



dt s 



2 



(15) 



with the invariant matrix elements squared and averaged (summed) over initial (final) colours 
and spins at LO in pQCD given by 

4 (t — ml) + (m — m^)^ + 2mls 



E 



qq-^cc 



I 9 ^2 



[ml — tj {ml — u) 8 (ml — tj {ml — u) — 2m^ 



{ml - i) {ml -u)- 2ml i^l + u) 2 ml{s- Ami) 



+ 



3 {ml - uf 3 (ml - Pj {ml - u) 

(ml — tj {ml — u) + ml (u — tj 

s (ml — t^^ 
(ml — t) {ml — u) + ml (i — u) 

6 ^ ^ .;^2_.2. ^ — ^- (16) 



s [m^ — 



Writing the four momentum of the incoming and outgoing particles as 



Pi = — (a;i,0,0,xi) 

P2 = — (a;2,0,0, -xa) 

Pc = {mTCOsh{yc),PT,0,mTsmh{yc)) 

Pc = (mTCOsh(?/g), -pr,0,mTsinh(?/g)) , (17) 
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the Mandelstam variables appearing in Eqs. ^ are given by 



s = 2m^ (1 + cosh(Ay)) 

i = — (1 + exp(— Ay)) 

u = ml - rriT {I + exp{Ay)) (18) 
Ay = Vc- Ve- 
in our calculation we use the GRV-LO ^ parton distributions in protons and pions, and 
apply a global factor of 2 — 3 in order to account for NLO effects. For the fragmentation 
function we use the Peterson function. 



3.2 Charmed meson production by recombination 

In the scenario described in for tt~ proton collisions, the annihilation of a m quark from 
the proton and the u quark in the pion would liberate the d of the pion which in turn re- 
combines to form a D~ and certainly not a D~^. On a similar base a p° proton collision will 
favor the production of and D~ over and depending on the quantum state of the 
colliding p° at the interaction point (see Fig. 3). 

The production of leading mesons at low p-r was described by recombination of quarks 
long time ago [0. In recombination models one assumes that the outgoing hadron is pro- 
duced in the beam fragmentation region through the recombination of the maximum number 
of valence quarks and the minimum number of sea quarks of the incoming hadron. The in- 
variant inclusive xp distribution for leading mesons is given by 

2E da''^^ f^F dxi dx2 



f^p dxi dX2 _ / XT-,/ N 

/ F2{Xi,X2) R2{Xi,X2,Xf) (19) 

JO Xi X2 



y/sa dxp Jo xi X2 

where Xi, X2 are the momentum fractions and F2 {xi,X2) is the two-quark distribution func- 
tion of the incident hadron. R2{xi,X2,xf) is the two-quark recombination function. 
The two-quark distribution function is parametrized in terms of the single quark distribu- 
tions. For recombination of D~ , 

F2 {Xi,X2) = PFd,u;val (^l) F^-sea {X2) (1 - Xi - X2) , (20) 

with Fg {xi) = XiQ (xi). We use the GRV-LO parametrization for the single quark distribu- 
tions in Eq. |20|. It must be noted that since the GRV-LO distributions are functions of 
X and Q^, our F2 (xi, X2) also depends on Q^. 
The recombination function is given by 

R2 {xu,d, Xc) = a ^"'^^^ ^ (^Xu,d + Xc- xf) , (21) 

X F 
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with a fixed by the condition dxF{^/<j)da'^^'^/dxF = 1- 

Some time ago V. Barger et al. |Tl[] explained the spectrum enhancement at high xp in 
production assuming a hard momentum distribution of charm in the proton. Here we will 
also take a QCD evolved charm distribution, of the form proposed by V. Barger et al. [|ll] 



xc{x,{Q'^)) = Nx\l-xf, (22) 

with a normalization N fixed to 

j dx ■ xc{x) = 0.005 (23) 

and I = k = 1. With this values for / and k one tries to resemble the distribution of valence 
quarks. In contrast with the parton fusion calculation, in which the scale of the inter- 
action is fixed at the vertices of the appropriated Feynman diagrams, in recombination the 
value of the parameter should be used to give adequately the content of the recombining 
quarks in the initial hadron. We used = 4m^. 



4 and D^{D^) total production 

The inclusive production cross section of a D meson is then obtained by adding the contribu- 
tion of recombination, Eq. |19|, to the QCD processes from direct photon-gluon interaction, 
quark anti-quark annihilation and gluon-gluon fusion from the hadronic component of the 
photon, i. e. 



da'"\D-) fda^s fdal^M , da- 



dxp \ dxp \ dxp dx 



TCC 

+ a { 6^1^ + ^rrZ^ ) ) (24) 



da'"\D+) (da^^ , ^ , dayj,^, 



^rec 



+ a h^^^}^ + (ilirZi^ (25) 



with 



dxp \ dxp \ dxp dx 



dxp dxp dx 



+ -IZ- (26) 



and A^i = , N2 = -^^J^^^ the parameters a, b, c and d depend on the contribution of 

each process to the total cross section. They are fixed in such a way that the differential 
cross section is well described, before calculating the asymmetry. The resulting inclusive 
D production cross section d(j^°^/dxF (shown in Fig. 4), is used then to construct the 
asymmetry defined in Eq. |1|. 

The values obtained for the different contributions in Eqs. 24 and 25 are in reasonable 
agreement with what one would expect. The photon fluctuation to a vector meson is of the 
oder of 1 %. Approximately 96 % of the total cross section comes from the photon gluon 
process. The contribution due to recombination goes from about 1 % (for D~) to about 3 % 
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(for D^). Fig. 5 shows the model prediction for the D^, production asymmetry together 
with the experimental results from the E687 Collaboration [|1|. The two curves correspond 
to the conventional GRV function distribution in the resolved photon and to the distribution 
proposed by Barger, et al. where a hard charm component has been assumed. 



5 Conclusions 



In an earlier work |12| the production asymmetry of Ac was described using the same recom- 
bination scheme used here. In hadroproduction the presence of a diquark in the initial state, 
plays an important role in production. In photo production, however, the production 
mechanism is somewhat different and the asymmetry is much smaller. The parameters used 
here are in reasonable agreement with what is physically expected and with the values used 
in a previous study of production asymmetries in vr proton collisions [^]. Changing the values 
of these parameters may improve the description of the asymmetry but then they may loose 
meaning in the frame of the asymmetry obtained for hadro production. 

Other experimental results for the asymmetry defined as: 

^^^^^ = a{D^)^a{D-) ^'^^ 

are A = -0.0384 ±0.0096 in xp > 0.0 at photon energies of 200 GeV and A = -0.0196 ± 
0.0147 in xp > 0.2 at photon energies of 80-230 GeV, average energy 145 GeV [|T3[. The 
NA14 collaboration studied the asymmetry 

^^^^^ - + DO) + a{D- + D^) ^^^^ 

and obtained A = -0.03 ± 0.05 xp > 0.0 at photon energies of 40-140 GeV and 100 GeV in 

average 



All these results are in good agreement with each other but, the statistical errors are still 
large. 

In ref. ^ experimental data are compared to a model based on string fragmentation. 
This model gives a larger asymmetry than the one obtained in our approach. The descrip- 
tion obtained there is in much better agreement with the experimental results. However, 
more precise measurements are needed before one can draw a final answer. 

One would expect that with increasing energy the resolved photon component increases 
giving rise to a larger contribution of the recombination mechanism which in turn would 
produce a larger production asymmetry. The HERA experiments should therefore be able 
to see a production asymmetry. However larger energies of the photon means also smaller 
values of x for the quarks that participate in the interaction and the hard charm component 
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is expected to play a minor role at small x's. i e., the asymmetry should look rather as 
the one obtained from conventional densities. New experiments will have soon more precise 
measurements of the asymmetries. This new results will give us a better understanding of 
the underlying production mechanisms. 
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Figure Captions 



Fig. 1: According with perturbative QCD, photon-gluon fusion is the main process in charm 
photo-production. 

Fig. 2: The resolved photon may interact with the proton via QCD, i.e. quark-antiquark 
annihilation and gluon-gluon fusion as shown in (a,b). A non perturbative interaction 
(c) may favor charm over anti-charm mesons production. 

Fig. 3: After a fluctuation of the photon to a p° vector meson, the interaction with the proton 
may ocurr in one of the two states. In any case the valence quark in the p would be 
released once the antiquark annihilates with a quark of the proton. 



Fig. 4: Total cross section as a function of xp- Experimental results from ||T3[ and theoretical 



calculation as in Eqs. El and p5| using the GRV distributions. 



Fig. 5: Measured production asymmetry for D~ and from [1]. The curves show the model 
result in which a hard charm component (dashed line) and GRV-LO (solid line) in the 
pion has been considered. The horizontal line at A{xf) = is for reference only. 
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Figure 1: 
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Figure 3: 
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Figure 4: 
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